Microscopic calculation of neutrino mean free path inside hot neutron matter 
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We calculate the neutrino mean free path and the Equation of State of pure neutron matter at 
finite temperature within a selfconsistent scheme based on the Brueckner-Hartree-Fock approxima- 
tion. We employ the nucleon-nucleon part of the recent realistic baryon-baryon interaction (model 
NSC97e) constructed by the Nijmegen group. The temperatures considered range from 10 to 80 
MeV. We report on the calculation of the mean field, the residual interaction and the neutrino 
mean free path including short and long range correlations given by the Brueckner-Hartree-Fock 
plus Random Phase Approximation (BHF+RPA) framework. This is the first fully consistent cal- 
culation in hot neutron matter dedicated to neutrino mean free path. We compare systematically 
our results to those obtain with the DIP Gogny effective interaction, which is independent of the 
temperature. The main differences between the present calculation and those with nuclear effective 
interactions come from the RPA corrections to BHF (a factor of about 8) while the temperature 
lack of consistency accounts for a factor of about 2. 
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I. INTRODUCTION 

Neutrinos play a crucial role in the physics of super- 
nova explosions [1] and in the early evolution of their 
compact stellar remnants 0, 01 • During the collapse of 
the pre-supernova core, a large number of neutrinos is 
produced by electron capture process. The mean free 
path A of these neutrinos decreases as the radius of the 
newly formed neutron star shrinks from about 100 km to 
about 10 km, becoming smaller than the stellar radius 
when density reaches a critical value (neutrino trapping 
density) . Under these conditions neutrinos are trapped in 
the star. Neutrino trapping has a strong influence on the 
overall stiffness of the dense matter Equation of State 
(EoS) 0,01, being the physical conditions of the hot and 
lepton-rich newborn neutron star substantially different 
from those of the cold and deleptonized neutron star. 

The scattering of neutrinos on neutrons is mediated by 
the neutral current of the electroweak interaction. In the 
non-relativistic limit for neutrons, the mean free path of 
a neutrino with initial energy E v is given by @ 

\- x (Eu,T) = ^J dk 3 (4(l + cos0) S<°>(?,T) 

+4(3-cos0) 5W(g,T)^ , (1) 

where T is the temperature, Gf is the Fermi constant, cv 
(ca) the vector (axial) coupling constant, fci = (E v , ki) 
and &3 are the initial and final neutrino four-momenta, 
q = fci — &3 the transferred four-momentum, and cos 9 = 
ki ■ k 3 . In the following, we impose the average energy 
E v = 3T [3. The dynamical structure factors S^ s \q,T) 



describe the response of neutron matter to excitations 
induced by neutrinos, and they contain the relevant in- 
formation on the medium (cf Eq. Ilfijl . The vector (axial) 
part of the neutral current gives rise to density (spin- 
density) fluctuations, corresponding to the S = (S = 1) 
spin channel. 

In this paper, we report on calculations of the mean 
free path of neutrinos in pure neutron matter under 
various conditions of density and temperature. A mi- 
crocopic framework based on the Brueckner-Hartree- 
Fock (BHF) approximation of the Brueckner-Bethe- 
Goldstone (BBG) theory is employed to describe consis- 
tently both the EoS and the dynamical structure factors 
of neutron matter including finite temperature effects. 

The paper is organized in the following way. A brief 
review of the BHF approximation at zero temperature 
and its extension to the finite temperature case is pre- 
sented in Sec.|n| The Landau parameters Fq, F\, Go and 
G\ are calculated in Sec. IIIII Section IIVI is devoted to 
the calculation of the dynamical structure factors and the 
neutrino mean free path. Finally, the main conclusions 
of this work are drawn in Sec. El 



II. EOS AT FINITE TEMPERATURE 

Our many-body scheme is based on the BHF approx- 
imation of the BBG theory. It starts with the construc- 
tion of the neutron-neutron G-matrix, which describes in 
an effective way the interaction between two neutrons in 
the presence of a surrounding medium. It is formally ob- 
tained by solving the well known Bethe-Goldstone equa- 
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FIG. 1: Mean field for several temperatures as a function of 
the momentum k in units of fef. The result of Brueckner- 
Hartree-Fock calculation is represented with filled circles 
while the mean field given by the DIP Gogny effective in- 
teraction is represented with empty squares. 



tion, written schematically as 



G(lo) = V + V 



Q 



U-E1-E2 



a] 



(2) 



In the above expression V denotes the bare interaction, 
Q is the Pauli operator which allows only intermediate 
states compatible with the Pauli principle, and oj, the 
so-called starting energy, corresponds to the sum of non- 
relativistic single-particle energies of the interacting neu- 
trons. The single-particle energy E is given by 



E{k) = 



h 2 k 2 

2m 



U(k) , 



(3) 



where the single-particle potential U(k) represents the 
mean field "felt" by the neutron due to its interaction 
with the other neutrons of the medium. In the BHF 
approximation U(k) is given by 



U(k) = Re n(k') (kk' G{u> = E{k) + E{k')) 

k' 

where 



kk' 



n(k) 



l,if k< k F 
0, otherwise 



A 

(4) 
(5) 



is the corresponding occupation number and the matrix 
elements are properly antisymmetrized. The resulting 
single-particle potential is shown toghether with the cor- 
responding one for the DIP Gogny || effective interac- 
tion in Fig.Q]for p = 0.25 fm~ 3 (k F = 1.95 fm~ 3 ). We 
note here that the so-called continuous prescription has 
been adopted for the single-particle potential when solv- 
ing the Bethe-Goldstone equation. As shown by the au- 
thors of Refs. the contribution to the energy per 
particle from three-body clusters is diminished in this 
prescription. We note also that the present calculations 
have been carried out by using the nucleon-nucleon part 
of the recent realistic baryon-baryon interaction (model 
NSC97e) constructed by the Nijmegen group The 
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FIG. 2: Energy per particle, free energy per particle and en- 
tropy per particle as a function of the number density for 
several temperatures. 



total energy per particle, E/A, is easily calculated once 
a self-consistent solution of Eqs. 10-0 is achieved 



E 
A 



h 2 k 2 

2m 



-U(k) 



(6) 



The many-body problem at finite temperature has 
been considered by several authors within different ap- 
proches, such as the finite temperature Green's function 
method ^3], the thermo-field method or the Bloch- 
De Domicis (BD) diagrammatic expansion 0]. The lat- 
ter, developed soon after the Brueckner theory, repre- 
sents the "natural" extension to finite temperature of the 
BBG expansion, to which it leads in the zero temperature 
limit. Baldo and Ferreira ^3] showed that the dominant 
terms in the BD expansion were those that correspond 
to the zero temperature BBG diagrams, where the tem- 
perature is introduced only through the Fermi-Dirac dis- 
tribution 



f(k,T) 



1 + exp ([£(*;, T)-/i„(T)]/T) 



(7) 



Therefore, at the BHF level, finite temperature ef- 
fects can be introduced in a very good approximation 
just replacing in the Bethe-Goldstone equation: (i) the 
zero temperature Pauli operator Q = (1 — «i)(l — 112) 
by the corresponding finite temperature one Q(T) = 
(1 — /i)(l — /2), and (ii) the single-particle energies E(k) 
by the temperature dependent ones E(k 1 T) obtained 
from Eqs. II3I4J1 when n(k) is replaced by f{k,T). 

In this case, however, the self-consistent process im- 
plies that together with the Bethe-Goldstone equation 
and the single-particle potential, the chemical potential 
of the neutron, /x„(T), must be extracted at each step of 
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U(k=0) m */m 




FIG. 3: Quadratic approximation of the single-particle en- 
ergy: U(k = 0) is the constant mean field and m* is the 
effective mass. See Fig. U for the legend. 
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FIG. 4: Landau parameters deduced for T=10 MeV (dotted 
line) and 80 MeV (dashed line) using the Nijmegen potential. 
It is compared with the one deduced from the DIP Gogny 
interaction (empty squares). 



the iterative process from the normalization condition 



£/(*,t). 



(8) 



This is an implicit equation which can be solved nu- 
merically by e.g., the Newton-Rapson method. Note 
that, now, also the Bethe-Goldstone equation and single- 
particle potential depend implicitly on the chemical po- 
tential. 

Once a self-consistent solution is obtained the total free 
energy per particle is determined by 



A A A 



(9) 



where E/A is evaluated from Eq. © replacing n(k) by 
f(k,T) and the total entropy per particle, S/A, is calcu- 
lated through the expression 



r -i^[/(i,r)k(/(fc,r)) 

k 

+ (l-/(fc,T))ln(l-/(fc,T))] . 



(10) 



Results for the quantities E/A, F/A and S/A are shown 
in Fig.|2lfor several densities and temperatures. 

The k- momentum dependence of the single-particle en- 
ergy is usually approximated by a quadratic function 
which consist in a constant term, the mean field, and a 
squared momentum dependent term, related to the neu- 
tron effective mass to*, reading 



E eS (k) 



H 2 k 2 
2m* 



U{k = 0) 



where 



h 

TO* 



l dU(k ) , 
k 5k 



k=k f 



(11) 



(12) 



We finish this section by showing in Fig.Elthe mean field 
and the effective mass as a function of the density, for 
several temperatures. As the temperature increases, the 
value of the mean field also increases. One can remark 
that the effective mass is nearly independent of the tem- 
perature up to T=80 MeV. 



III. LANDAU PARAMETERS Fo, Fx, Go AND Gi 

The Landau parameters Fo, Go are related to the in- 
compressibility modulus K and the magnetic susceptibil- 
ity x respectivly through 



K = 9p^ = ^ + F ° 



dp 2 



Nn 



X 



1 d 2 F/V 1 + G 



(up) 2 dp 2 s 



XF 



(13) 



(14) 



where F/V = (F/A)p is the free energy density, No = 
m*kp /ir 2 h 2 is the density of states, \f is the magnetic 
susceptibility of a free Fermi gas and p s = p-\ — pi is 
the spin asymmetry density (i.e., the difference in the 
densities of neutrons with spin up and spin down) . The 
parameter F\ is deduced from the effective mass accord- 
ing to the relation 



m 



1 + f 

to 3 



(15) 



There is no simple relation for the Landau parameter Gi 
that can be deduced from thermodynamical properties 
or general relations like the forward scattering sum rule. 
Gogny effective interaction predicts that G\ is close to 
zero at all densities (cf Fig. HJ, hence, we will assume in 
the following that Gi=0. 
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Finally, the Landau parameters we obtain in pure neu- 
tron matter are shown on Fig. Q] The absolute values of 
F and Go increase with the temperature while F\ is 
nearly constant. The main differences between the Lan- 
dau parameters calculated with the BHF approximation 
and those calculated with the Gogny effective interaction, 
are present in the spin fluctuation channel. In the latter 
case, for high densities, the Landau parameter Go falls to 
— 1 because of a ferromagnetic instability while Go 
increase with density for the BHF calculation, since no 
such an instability is found in this case [ItI ]. 



IV. DYNAMICAL STRUCTURE FACTORS 
AND NEUTRINO MEAN FREE PATH 

We obtain the dynamical structure factors from 
the imaginary part of the response function \^ (l, T) in 
pure neutron matter as fli| 



0,005 



S (5) (q,qo,T) 



7T 1 - exp(-q /T) 



3m X ( S )(q,q ,T). 



(16) 

In the Brueckner-Hartree-Fock (BHF) approximation 
the response function is given by 



3mv (S) fa a T) — ^ In 
X BHF (q, q , J- ) - — 111 



1 



3 (A+q /2)/T 



1 + e (A-q /2)/T 

'(17) 

This expression 



where A = \i n - m*/2(q /q) 2 - q 2 /& 
is valid for both positive and negative transferred energies 
qo. The real and imaginary part of the response function 
have been obtain within the quadratic approximation for 
the single-particle energy (see Eq. dJ). The influence of 
the in-medium interaction appears through the effective 
mass to* and the chemical potential fi n , which introduce, 
besides a density dependence, an additional temperature 
dependence in the response function. 

We have also included the long range correlations 
within the BHF+RPA scheme. The particle-hole inter- 
action is approached by a Landau form, containing only 
I = multipole (RPA I = 0), or I = 0, 1 multipoles (RPA 
I = 0, 1) [13 • We stress that in all cases the particle- 
hole interaction has been consistently obtained from the 
particle-particle interaction used to describe the EoS. 

The dynamical structure factors, within the BHF and 
BHF+RPA scheme, are displayed in Fig. and Fig. as 
a function of the transferred energy qo and for p=0.25 
fm~ 3 , T=10 MeV and q = 10 MeV. In Fig. we show 
the dynamical structure factors for the density (S = 0) 
and spin-density (S = 1) channels. We compare the re- 
sponse function obtained with the Brueckner input (filled 
circle) to the result with the DIP Gogny effective inter- 
action (empty square). The two columns correspond to 
two truncations in the BHF+RPA calculation which con- 
sist to include only the I = Landau parameters or the 
I — 0, 1 Landau parameters. As expected from the values 
of the Landau parameters (cf Fig.^J, the main differences 
between the approaches are present in the spin channel 




5 
q [MeV] 



FIG. 5: Dynamical structure factors at the BHF approxima- 
tion as a function of the transfered energy qo, and for p = 0.25 
fm, T=10 MeV and q = 10 MeV. 



RPA 1=0 



RPA 1=0,1 




q„ [MeV] 



q n [MeV] 



FIG. 6: 
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Dynamical structure factors, for the channel S = 
1, as a function of the transfered energy qo, and for 



fm 



T=10 MeV and q=10 MeV. 



(S = 1). The spin sound is weakly present for the DIP 
Gogny effective interaction while it seems very impor- 
tant for the Brueckner calculation. The pronounced zero 
sound in the spin-density channel at zero temperature is 
damped when the temperature increases, and these dif- 
ferences disappear at high temperature. 

The density and temperature dependence of the neu- 
trino mean free path is shown in Fig. We have also 
performed a comparison with Gogny prediction. The be- 
haviour of the mean free path follows approximatively 
the expected law in T~ 2 0|. At low temperature, the 
discrepancies between the two interactions are essentially 
due to the density dependance of the effective mass (cf 
Ea. ll7p . Moreover, as the effective mass is independent of 
the temperature, a calculation of the in-medium interac- 
tion at finite temperature is not necessary. Hence, within 
BHF scheme, the mean free path is sensible to the tem- 
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FIG. 7: The neutrino mean free path deduced from BHF 
response function as a function of the density for T=10 MeV 
(left) and T=80 MeV (right). The energy of the incoming 
neutrino is E v = 3T. 
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FIG. 8: Effects of BHF+RPA I = 0, 1 for Brueckner (filled 
circle) and Gogny interactions (empty square). The solid line 
stands for a consistent calculation at T=10 MeV, the dashed 
line is for T=80 MeV but with the EoS at T=10 MeV (con- 
sistency broken) and the dotted line is for T=80 MeV (con- 
sistent). The energy of the incoming neutrino is E v = 3T. 



perature essentially through the explicit T-dependance of 
the response function and the effect due to the finite tem- 
perature EoS is negligable. At high temperature, the 2 
interactions gives very similar results because the effects 
of the medium are reduced (E v = 3T). 

The effects of the residual interaction, i.e., of RPA cor- 
relations, can be seen in Fig. |H1 where we have repre- 
sented the ratio R = X BliF+RPAl=0 - 1 / A BHF of the neu- 
trino mean free path within BHF+RPA I = 0, 1 to BHF 
as a function of the density. The solid line stands for 
T=10 MeV and the dotted line for T=80 MeV. The ra- 
tio R with Nijmegen potential is about 8 at high den- 
sity while it is only 1 with Gogny interaction. This ef- 
fect have also been shown for zero temperature EoS with 
three-body forces j^J. This is due to the presence of spin 
instabilities at high density in the case of Gogny interac- 
tion [3 EH. These instabilities lead to the divergence of 
the dynamical form factor ^| > hence the neutrino mean 
free path goes to zero. This illustrate the sensibility of 
the neutrino mean free path with the onset of instabili- 
ties. Here, it reduces the mean free path by a factor 8 at 
high density. 

In order to understand the interplay between the ex- 
plicit T-dependance of the response function (cf Eq. I17|) 
and the Landau parameters calculated at finite temper- 
ature, we have broken the consistency of the approach. 
In Fig. the solid line stands for a consistent calcu- 
lation at T=10 MeV while the dashed line stands for 
the calculation of the neutrino mean free path using low 
temperature Landau parameters (T=10 MeV) but cal- 
culating the dynamical response function at T=80 MeV. 
Hence, by comparing the solid and the dashed line cal- 
culated with the same particle-hole residual interaction, 
we see that the increase of temperature reduces the ef- 
fect of the correlations. As an illustration, the increase 
of the temperature suppress the collective modes. The 



increase of temperature induce also the increase of the 
incoming neutrino energy (E u = 3T). The lack of con- 
sistency between finite temperature EoS and the dynam- 
ical response function leads to an underestimation of the 
neutrino mean free path of about 30-50%. The effect of 
restauring the temperature consistency for T=80 MeV 
is illustrated by comparing the dashed line and the dot- 
ted line (consistent calculation for T=80 MeV). As it has 
been shown on Fig. 0] the Landau parameters Fq and 
Go increase with the temperature. The correlations be- 
come more important and the ratio R increase. Finally, 
the temperature induce two opposite effects: it increases 
the particle-hole interaction, but it decreases the effects 
of these correlations in the calculation of the dynamical 
structure factors. These two effects tend to compensate 
each other and the ratio R, at low temperature, is close 
to the ratio at high temperature. 



V. CONCLUSIONS 

The purpose of this article is to study consistently the 
effects of the temperature to compute the neutrino mean 
free path in dense neutron matter. Both the EoS and 
the response function have been computed, including the 
BHF+RPA correlations. We have shown that the effec- 
tive mass is nearly independent of the temperature. 

We have shown the interplay between the effects of 
temperature for the calculation of the EoS (the mean 
field and residual interaction deduced from it) and the 
explicit temperature dependance of the dynamical struc- 
ture factor. The first one increases the correlations with 
the temperature while the second one decreases it. Fi- 
nally, these two opposite effects approximatively compen- 
sate each other and the ratio R is nearly constant for the 
range of temperature 10-80 MeV. 
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The temperature consistency does not change the gen- 
eral behaviour of the neutrino mean free path. The lack 
of consistency leads to an underestimation of the neutrino 
mean free path of about 30-50% for T=80 MeV. 
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